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Blazar Flares from Compton Dragged Shells 

Omri Golan and Amir Levinson ^ 

ABSTRACT 

We compute the dynamics and emission of dissipative shells that are subject to a 
strong Compton drag, under simplifying assumptions about the dissipation mechanism. 
We show that under conditions prevailing in blazars, substantial deceleration is antici¬ 
pated on sub-parsec and parsec scales in cases of rapid dissipation. Such episodes may 
be the origin of some of the flaring activity occasionally observed in gamma ray blazars. 
The shape of the light curves thereby produced reflects the geometry of the emitting 
surface if the deceleration is very rapid, or the dynamics of the shell if the deceleration 
is delayed, or initially more gradual, owing, e.g., to continuous injection of energy and 
momentum. 

Subject headings: galaxies: active - quasars: general - radiation mechanism: nonthermal 
- X-rays:galaxies 


1. Introduction 

The broadband spectrum observed in blazars is dominated by beamed emission produced in 
relativistic jets that emanate from the central black hole. These jets propagate through a dense 
radiative environment, and interact with seed photons that are supplied by extended radiation 
sources, notably the accretion disk around the black hole, gaseous clouds in the broad line region 
(BLR), and a dusty molecular torus located at larger scales (see, e.g., Joshi et al. 2014 for a recent 
account). This interaction should affect the dynamics and emission of the jet. The primary concern 
of most previous works (e.g., Dermer & Schlickeiser 1993, Sikora et al. 1994, Blandford &: Levinson 
1995, Ghisellini & Madau 1996; see also Levinson 2006 and references therein) has been the effect 
of this interaction on the observed spectrum, and simple emission models that ignore dynamical 
effects and complex structures have been constructed for this purpose, although some recent works 
incorporate more realistic models for the dynamics of the emitting plasma (e.g., Joshi &: Bottcher 
2011, Joshi et al. 2014, and references therein). In ERG models, the high-energy component of 
the spectral energy distribution (SED) is attributed to inverse Gompton scattering of ambient seed 
photons by non-thermal electrons accelerated in dissipative regions inside the jet. This process, 
unlike synchrotron emission, gives rise to a loss of linear momentum of emitting fluid elements 
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and a consequent radiative drag that tends to decelerate the bulk flow. In certain circumstances, 
explored below, this can lead to rapid, large amplitude flares. 

Rapid variability over the entire electromagnetic spectrum is a characteristic property of 
blazars. Episodic gamma-ray emission with flares durations of hours to weeks is quite typical 
to many gamma-ray blazars, with the most extreme activity recorded in the sub-class of TeV 
AGNs, e.g., Mrk 421, Mrk 501, PKS 2155-304. The observed variability time imposes a stringent 
constraint on the maximum size of the emission region, Ar, which in extreme cases is inferred to 
be of the order of the gravitational radius of the putative black hole. The naive expectation is that 
the large amplitude short duration flares seen in gamma-ray blazars originate from small radii, as 
the fraction of jet energy that can be tapped for production of 7 rays scales as 7 ~ {ArfOjrem)^, 
where Oj is opening angle of the jet and rem the emission radius, and is small for rem » Ar. 
Some of the variable gamma-ray flux may be attributed to a sparking gap at the base of the jet, 
as proposed for M87 (Levinson 2000, Neronov &: Aharonian 2007, Levinson &: Rieger 2011), how¬ 
ever, in powerful blazars the spectrum emitted from the gap is unlikely to extend beyond a few 
GeV, owing to the large pair-production opacity. This, and the indications of correlated emission 
at much lower energies (radio-to-x rays) strongly suggest that the jets are important sources of 
episodic gamma-ray emission. It has been argued that in some cases the observations favor models 
in which the variable gamma-ray emission originates from small regions locate at large radii (e.g., 
Sikora et al. 2008, Agudo et al. 2011). Such events may be produced by a converging shock in a 
reconfinement nozzle (Bromberg &: Levinson 2009) or magnetic reconnection in minijets (Giannois 
2013). However, recent analysis (Nalewajko et al. 2014) challenges the far dissipation scenarios, 
indicating typical emission radii in the range 0 . 1-1 pc. As shown below, on these scales the radiative 
drag can be substantial. 

In this paper we consider the effect of Gompton drag on the dynamics and emission of dissipa¬ 
tive fluid shells. We show that if dissipation of the bulk energy commences not too far out, the blob 
experiences strong deceleration that leads to large variation of the observed flux emitted from the 
blob by virtue of the change in its bulk Lorentz factor. The effect of Gompton drag on the dynamics 
of a relativistic jet has been considered earlier by several authors under different assumptions (e.g., 
Phinney 1987, Li et al. 1992, Sikora et al. 1996). Our method is similar to that presented in 
Sikora et al. (1996), however, we focus on short events that may lead to rapid flares, and compute 
the Lorentz factor profiles and the resulting gamma-ray lightcurves for a range of conditions. For 
typical ambient luminosities, the duration of flares produced by this mechanism is on the order of 
the characteristic size of the emitting blob. Thus, flare durations as short as the dynamical time 
of the central engine can naturally be accounted for in this model, provided that the bulk energy 
can be dissipated at a large enough rate. A preliminary account of the decelerating shell model is 
given in Levinson (2007). Here, we present an elaborated analysis of the dynamics of the flow, and 
also compute the resulting light curves. The construction of this model was originally motivated 
by the apparent discrepancy inferred in TeV blazars between the relatively large Doppler factors, 
(Id ~ 30 — 50, required to avoid strong attenuation of the VHE flux emitted from the inner regions. 
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and the much lower values, ~ a few, inferred from superluminal motions and source statistics 
(Georganopoulos &: Kazanas 2003, Levinson 2007). However, the question of how the Compton 
drag affects the dynamics and emission of dissipative outflows is of general interest, and is relevant 
essentially to all blazars. 


2. The model 

In the simplified model invoked here, a blob is ejected from a central engine of size r* = 2GM/(? 
and accelerated to a Lorentz factor Fq. When it reaches some radius ~ F^r^ ~ IO^^FIqMq cm, 
dissipation suddenly starts, e.g., via formation of shocks or explosive conversion of magnetic energy, 
leading to continuous particle acceleration. Inverse Compton scattering of ambient photons then 
leads to a strong radiative drag that tends to decelerate the blob. In certain circumstances, this 
drag may be compensated by injection of energy and momentum into the blob by some external 
agent, delaying the deceleration of the emitting plasma. The deceleration may be delayed also in 
cases where the dissipation commences very close to the central engine, where the radiation field is 
dominated by direct illumination from the disk and is highly unisotropic, as in this zone the blob 
may propagate at the equilibrium Lorentz factor until reaching scales where the external radiation 
field is roughly isotropic (e.g., Sikora et al. 1996, Vuillaume et al. 2014). In the present model we 
shall refer to this case as delayed deceleration. In reality, the structure and velocity of the emitting 
plasma are expected to be non-uniform by virtue of confined dissipation and rapid cooling. In what 
follows we ignore such complications and compute the dynamics and emission of the blob under 
simplifying assumptions about the microphysics of dissipation. 

The energy and momentum fluxes of the emitting fluid are given by, 

= (u; + HV47r)F2/3r = u^T^/Ir, (1) 

= {w + HV 47 r)F 2/32 + ^ (2) 

where w = p + e + p \s the proper specific enthalpy, p, p and e are the pressure, proper density 

and proper internal energy, respectively, B is the proper magnetic field, and F = (1 — is 

the bulk Lorentz factor. In terms of the total power, Pj = , where 9j is the half opening 

angle of the flow, we have 

^ = ^g2^2cr2 = O.lP,-44(0jrri7)-2 erg cm-3, (3) 

for p-p = 1, r = 10 ^^ri 7 cm and Pj = 10 ^^Pj 44 erg s“^. 

For simplicity, we assume that the ambient radiation field is roughly isotropic in the frame 
of the central engine, and has a luminosity Lg. This quasi-isotropic radiation field is contributed 
by scattering and reprocessing of the central UV radiation by gas in the broad line region, and 
by emission from a dusty torus (e.g., Joshi et al. 2014). Here, Lg represents the sum of these 
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components, and is a fraction of the total luminosity of the continuum source. The radial profile 
of the ambient intensity intercepted by the jet depends on the geometries of the broad line region 
and the dusty torus (e.g., Joshi et ah, 2014). Denoting x = r/r^, we express it as 


Is{es■,^i,r) = Ks{r)F{es)., Hs{r) 


L, 


lOvr^r^c 


fsix), 


(4) 


with the normalization F{€s)d€s = 1, such that the total energy density at radius r is Us{r) = 
4:TrKs{r). For central emission, fs{x) = x~‘^ = (r^/r)^. For a flat profile within some scale ri > r^, 
fsix) = {vd/ri)'^ at X < ri/vd- For the spectrum of the ambient radiation field we adopt 


F(es) oc < 


^1/2 

ts 

- 3/2 

es 


0 

\ 


10 ® < Cs/rUeC^ < 10 
10 -“^ < €s/meC^ < 10 -^ 
else 


that mimics a typical soft spectrum. 


(5) 


In the rest frame of the blob the intensity is given by 

and the comoving energy density by 


u' = 2tt r, 


.de',d/i' = 2TTKsir) J ^ ^ + /3r/3). 


( 6 ) 

(7) 


From Equation ([6]) it is seen that the comoving intensity peaks sharply around the direction opposite 
to the blob velocity, viz., fi' = —1. To simplify our calculations we approximate the comoving 
intensity as a beam moving in the direction /i' = —1: 


= dsir)Fie'Jil + I3rfj.'))5il + fj.'), ( 8 ) 

with rjsir) = Ks(r)2r(l + /3r)“^(l + /dr/3) — 4FKs(r)/3. It can be readily verified that with this 
choice the comoving energy density satishes Equation (|7]). 

We further suppose that the electron distribution function is isotropic in the fluid rest frame 
and can be approximated as a power law: = Keir)'^~^\ 7i < 7 < 72 , where me(?^ is 

the corresponding electron energy, as measured in the comoving frame. We define f,eir) to be the 
fraction of the total internal energy carried by the relativistic electrons, that is, = n/, where 
Ug = f meC^jdng, subject to the boundary condition ^e(r < rd) = 0. In terms of this parameter 
Ke = ^e(ujjmeC^)(2 - q)/(j 2 ~'^ “ 7 i~'*')> and we note that ( 72 “'' - j^~'^)/(2 - q) = ln( 72 / 7 i) in 
the limit g —>■ 2. We suppose that reaches a maximum value over some characteristic length 
scale that depends on the microphysics of the specific dissipation mechanism. This length scale is 
typically of the order of the gyroradius of accelerated electrons, which is shorter than the cooling 
time of the highest energy electrons. Thus, we take to be constant at r > r^. 
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The maximum Lorentz factor of the electron distribution, 72 , is likely to be limited by cooling. 
In the limit of Bohm diffusion the acceleration time is = {'yniec)/{r]acc^B). Equating with the 
Compton cooling time, = jmeCp/Pcom = and adopting ijacc = 0.1, we estimate 

r72 ^ 3 X 

in terms of the fraction = 0.1 ^b_i = u'^ = Here Bq = B{r = r^) is the magnetic 

field at the onset of dissipation. 


2.1. Dynamics 


We consider only situations in which dissipation is rapid enough to produce high enthalpy 
inside the blob on a timescale comparable to the radiative time. Then, one can adopt a simple 
prescription in which the dissipation is formally treated as an initial condition (at r = r^) of the 
flow equations. The magnetization of the flow is then expected to be low, so that the effect of 
the magnetic field on the dynamics of the system can be neglected. This may not apply to cases 
were the dissipation time is much longer than the light crossing time of the shell. We suppose that 
the dissipated energy is redistributed in a way that a fraction of the total dissipation energy 
is injected in the form of a power law electron distribution, as explained above. Under these 
assumptions, the dynamics of the blob subsequent to the onset of dissipation is governed by the 
equation 

= s,", ( 10 ) 

where the source terms account for the radiative drag acting on the blob, and are given explicitly 
in Equations (jA15p - (lAlOp . Eor a conical expansion, the 0 and r components of Equation (IlOp can 
then be written, to order 0(r“^), in the form (see appendix [B]) 


d 1 r dp 

w— Ini + — 
dt dt 


- r-25 = S: - /3r50, 


9C 


^ d 2 d dp <9 , T. ci 


dt 


dr dC dC 


( 11 ) 

( 12 ) 


in terms of the coordinates t and Ci^p) = ^ 2 (i) — f, where d/dt = dt + Prdr is the Lagrangian 
derivative, and r 2 (t) = ct + 0(r“^) denotes the trajectory of the blob’s front. The interior of the 
blob encompasses the range 0 < C < A- 


Now, for the relativistic electron population invoked above < 7^/3^ >=< 7 


< Y >= 


1 


n' 


2dn'e, , 


r7max 

h^in d-f 


dj, 


1 , where 

(13) 


and to order 0(r ^) the source terms. Equations (lAlSp - ()A16p . are given by 


5 ° = --T^ < 72 > UsarY., 


( 14 ) 
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and 


5,^ = /3r5° + 5°/3r2 


(15) 


2.1.1. Uniform blob 

We consider first a uniform blob ( 9 ^ = 0 ) having a fixed length A in the frame of the central 
engine. Its motion is then characterized by a single Lorentz factor at all times, r((^,t) = r(t). We 
suppose that microphysical processes redistribute the dissipation energy in a way that a fraction 
of the total internal energy is uniformly injected in the form of a power law electron distribution 
inside the blob. For clarity, we ignore here the contribution of the thermal population to the 
Compton drag. As shown below, this is justified when f,eXl2fne<^ is larger than the thermal energy, 
here 72 is the upper cutoff of the electron distribution given in equation Q, and x is defined below. 
The inclusion of the thermal electrons in the source terms does not alter our results significantly, 
and at any rate, if the dissipation produces a population of relativistically hot electrons with a 
thermal energy rather than a power law distribution, it can be readily accounted for by 

taking x?e = 1 , 72 = 7t in Equation (fl^ below. 

We find it convenient to use the parametrization < 7^ > / < 7 >= X72) where < 7 > rUeC^ = 
tig/rig is the average energy of nonthermal electrons. For the power law energy distribution invoked 
above we have 

_ (2-g) [1 - (72/71 71 

^ (3-g)[l-( 72 / 71 ) 2 -'/] 72- ^ ^ 

With 72 >> 7i and g < 2 it gives x — ( 2 —g)/(3—g). For q = 2we have x = [ln(72/7i)]~^ > 0 . 1 , and 
for q = 2 . 5 , x — 1/71 /72• A reasonable choice for the minimum electron energy, adopted henceforth, 
is 71 = mp/nie. We further define the fiducial coordinate x = ct/rd, and write Ug = Usofs{x)/fso, 
denoting fgQ = fs{x = 1). Equations (fTTD and (fT^ with d(^p = 9 ^F = 0 can then be re-expressed as 


dx 

din 


ln(r X w) = -6a(x/xo)(72/72o)(r/ro)(/s//so)(4e/3'u;), 


wx 


— —dlnp = 0, 
w 


in terms of the constant 


a = 


3mpc2 


— 3 X J sO 


l/2^-l/2 A/2 

J S 


(17) 

(18) 

(19) 


The solution of these coupled equations depend on the equation of state and the assumptions about 
72 . Approximate analytic solutions can be obtained in the case q < 2 for which x — Xo is a good 
approximation. For a relativistically hot blob, we adopt the equation of state e = 3p. We then 
obtain 


^ln(r/x) = -a(72/72o)(r/ro)(/s//so)- 


( 20 ) 
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As a first example, let us take the maximum electron energy to be constant during the deceleration 
phase, that is., 72 = 720- The solution then reads: 


r(x) 


Tqx 

[l + agi{x)y 


( 21 ) 


where gi{x) = jy[x'fs{x')/fso]dx', with gi{x) = Inx for and gi{x) = x ^/2 for a fiat 

profile, fs/fso = 1- As a second example, we suppose that 72 is given by Equation ([9]), and that 
the magnetic field evolution is dictated by the ideal MHD limit, viz., B/Bq = ro/(xr). We then 
obtain the solution 

T{x) = ro[l - ag 2 {x)y, (22) 

with g 2 {x) = jy x'-^[fs{x')/fsoy/‘^dx'. It is worth noting that these analytic solutions hold only 
for times x at which the blob is relativistically hot, and r(x) >> 1. 

Figured] exhibits exact numerical solutions of Equations (fT71l - (fT8]) for poc^/po = 10 “^, 72 given 
by Equation Q with B/Bq = ro/(xr), intensity profile /* = and different choices of a. As 
a check, we obtained solutions also for different intensity profiles and different assumptions about 
72 , and found little differences in the Lorentz factor profiles for a given a in the regime of rapid 
deceleration. The key parameter that determines the dynamics is a. 


2.1.2. Internal shocks 


The model outlined above assumes uniform acceleration of electrons to nonthermal energies 
at all times. This assumption may hold in certain situations but not in general. In shocks, for 
instance, particle acceleration is confined to a region around the shock front of a characteristic size 
comparable to the gyroradius of the accelerated electrons, as measured in the shock frame. For the 
highest energy electrons it is roughly equal to the cooling distance (if not limited by escape). For 
lower energy electrons it may be even smaller. As the downstream fluid moves away from the shock, 
electrons of energy yrUeC^ cool over time = mecj{'^aTu'g), where tt' = dT^tt^/S is the comoving 
energy density of the external radiation (see Eq. CD. The distance traversed by these electrons in 
the black hole frame before cooling down is Ic = rt(.c = 3meC^/(4r7crrUs). Using Equation (fT^ 
then yields alc/rd = x^e << 1 for << 1- Consequently, IC scattering off nonthermal electrons 
should not lead to significant deceleration of the blob, unless the entire shock energy is converted to 
non-thermal electrons {x^e = !)• However, about half the shock energy is carried by a population 
of thermal electrons having an average comoving energy ^just downstream of 

the shock, where denotes the shock Lorentz factor (that is, the Lorentz factor of the upstream 
fluid measured in the shock frame). These electrons cool as they propagate away from the shock. 
The Lagrangian rate of change of the thermal energy in the downstream flow is governed by the 
equation 


djT 

dt 







SmeC 


( 23 ) 





To compute the shock structure one needs to solve Equations (fTT]) - (fT^ coupled to Equation ([231) . 
Such treatment is beyond the scope of this paper. To illustrate the effect of radiative drag on the 
shock we invoke the uniform shell approximation, that is, keep only the Lagrangian derivatives in 
Equations (fTT]) - (fT^ . Then, adopting < 7 ^ > meC^n'^/w = and w = Ap we obtain the rate of 
change of the bulk Lorentz factor: 


^in(r/x) = -l^r7T. 

dt drUeC 

With Us = Usofs{x)/fso, the solution of the coupled Equations (1231) and (l2^ reads: 


(24) 


r(x) = To 


X 


[1 + argrix)]'^/'^' 


(25) 


here gxix) = x'^^/'^lfsix')/fso]dx', and 


CXt 


ISrdaT^ojToUsO 

OrUeC^ 


— ^^ 0 ^ sO' 


(26) 


Since ToT^/j >> 1 signihcant deceleration is expected on sub-parsec scales. For instance, assuming 
a flat intensity prohle below the radius r; = 10 ^® cm, with ttgo = 10 “^ erg cm~^ at r < n, we 
obtain ot = 0.2rors/irrfi7. Thus, colliding shells having a Lorentz factor Tq > 10 will experience 
substantial deceleration. 


Now, the deceleration of the downstream plasma should lead to a gradual strengthening of 
the reverse shock and a weakening of the forward shock. Our preliminary calculations indicate 
substantial over-compression of the reverse shock already at ~ a few. In the frame of the central 
engine this translates to a deceleration of the entire shocked shell, as described qualitatively by the 
simple blob model outlined in the preceding section. The strengthening of the reverse shock should 
lead to enhanced dissipation rate, whereby the bulk energy of the unshocked shell is ultimately 
radiated away with high efficiency. Consequently, as long as the beaming cone of emission is 
narrower than the angular extent of the shell, the total flux observed will remain roughly constant. 
Thus, the delayed-deceleration model is relevant for the evolution of the total flux. However, the 
change in the Doppler factor resulting from the deceleration of the shocked shell leads to a change 
in the observed energy of scattered photons and this, in turn, can significantly alter the evolution 
of the SED. In particular, we anticipate different durations and times of peak emission of flares 
observed in different energy bands. The non-uniformity of the emitting plasma downstream of the 
shock adds complexity. 

A comprehensive analysis of emission from internal shocks in blazars is given in Joshi & 
Bottcher (2011) and Joshi et al. (2014) ignoring the effect of Compton drag on the dynamics of 
the shock. The decay of the emission in their model is due to cooling of the emitting electrons 
following shock crossing. This situation is well represented by the delayed deceleration model we 
adopt below. The neglect of Compton drag is justified only at radii where a-r < 1 in Equation ([261) . 
As explained above, the inclusion of Compton drag and non-uniform particle acceleration should 
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have a profound effect on the light curves. A complete treatment of the dynamics and emission of 
internal shocks that are subject to a strong Compton drag will be presented in a future publication. 


2.2. Inverse Compton emission 

As will be shown below, the shape of the lightcurves reflects the geometry and dynamics of 
the emitting material. In particular, time retardation associated with the curvature of the emitting 
surface sets a limit on the rise and decay times of the flare. Furtheremore, as mentioned above 
certain situations can be described by delayed deceleration of the blob, that can have an important 
effect. To illustrate such effects, we shall consider also cases in which the Lorentz factor remains 
constant at its initial value Fq up to some radius < r^ec < + ^/(l ~ /^r), and only then 

deceleration commences. In reality, the structure of the emitting zone is expected to be non- 
uniform in those cases, depending on the specific model. Those details may affect the resulting 
emission. Our purpose here is merely to illustrate how dynamical effects are imprinted in the 
lighcurves. For this purpose our simple treatment of delayed deceleration is sufficient. Moreover, 
this prescription also describes situations in which the radiative drag is too small to affect the 
dynamics of the shell, and the decay of the emission at the end of the dissipation phase is due to 
cooling of the emitting electrons. 

Since we are mainly interested here in the high-energy emission from the blob, we consider 
only the contribution of Inverse Compton scattering. The spectral evolution measured by a distant 
observer is computed as follows: For a given choice of parameters we hrst solve Equations (II7p - (II8p 
numerically to obtain the Lorentz factor profile F(r) in the frame of the black hole. Adopting the 
beam approximation, the intensity of the background radiation at any radius r is then transformed 
into the rest frame of blob using F(r) in Equation ([8]). The comoving intensity thereby obtained 
is used to compute the comoving emissivity of the scattered radiation, /r', r). The emissivity 

in the black hole frame is given by jsc{e,fi,r) = [r(l — /3rjsd^', fJ-' , r) , where fj, = {fJ-' + 
Pr)/{1 + Pr fJ-') and e = T{1 + ■ The scattered intensity emitted by the blob is obtained upon 

integrating the emissivity across the blob, taking into account the time delay between emission of 
photons from different locations: 

pAr 

Isc{€,fi,r)= jsc{€,ti,r - y)dy (27) 

Jo 

where the distance Ar is related to the blob’s length A through 

A = / [1-/3r{r')fAdr', (28) 

J r—Ar 

and is a function of r. In deriving Equation (j27|) we assumed that the emissivity is uniform inside 
the blob and vanishes outside it. Note that 

tobir) = f [1 - f3r{r')y]dr'/c 
Jrd 


( 29 ) 
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is the time measured by a distant observer viewing the blob at an angle 9 = arccos(^) relative to 
its direction of motion, so that formally A/c = tob{r) — toh{f — Ar). 

High energy photons emitted by the blob will be attenuated by pair production on background 
photons. Here we model this attenuation by an exponential cutoff at the corresponding pair pro¬ 
duction optical depth T^^{e,r). The latter is computed using Equations (3.1)-(3.3) in Blandford & 
Levinson (1995), and plotted in figure [2l The spectral flux measured by a distant observer viewing 
the blob at time tabs at an angle 9 satisfies 

oc / dy f jsc(e,/i, r - (30) 

Jo ^S(r—y) 

where the radius r is computed at the observed time tob using Equation ([29]1 . that is r = r{tob), 
and the inner integration is over the emitting surface of the blob at the retarded time (r — y)jc. 


In the immediate deceleration case, the rise and decay times of the observed flux are dominated 
by temporal delays associated with the curvature of the emitting surface. Eor a conically expanding 
shell the emitting surface is spherical, and Equation (j3n|] reduces to: 


jro 


ryo ry-2ir,y) 

{e,iy,tob)oc. dy jsc{e,f^,r-y){r-yfe~'^'^'^^^’''~y^y~^dy, 

Jo J/ii{r,y) 

H 2 {r, y) = max | cos 9j, 1-/ [/3f ^(y') - l]dy' \ , 

t r-y } 


with 


^i(r,y) = min < 1,1 - 


{1- H 2 ){r-y) - A 
r — y — A 


(31) 


(32) 

(33) 


where 9j is the opening angle of the flow, and yo is determined from the condition /ri(r, yo) = cos 9j. 

The evolution of the SED is shown in figure [3] in the case of immediate deceleration (left panel) 
and delayed deceleration with -|- A/(l — /3r) (right panel). The corresponding lightcurves 

are displayed in figure 0] at a photon energy e = 1 GeV. The right pannels exhibit the emissivity, 
and the left pannels the observed flux computed for these emissivities using Equation (1311) . The 
effect of the jet opening angle on the shape of the lughtcurve is more prominent in the non-delayed 
case, as seen in the upper left panel of figured! in the delayed case this effect is essentially negligible. 
These lightcurves are rather typical in the regime where the deceleration time of the blob, Vd/ica), 
is not much larger than its light crossing time A/(l — /3r). In this regime, the duration of the flare 
measured by a distant observer viewing the source at an angle smaller than the opening angle of 
the flow and the overall shape of the lightcurve depend on details, as seen in figured! 


In the case of immediate deceleration the lightcurve is asymmetric, with the rise and decay 
times determined by the curvature of the emitting surface. For a sufficiently large Lorentz factor, 
such that the beaming angle is smaller than the opening angle of the flow (i.e., Tg^j > 1), these 
times, as measured by a distant observer, satisfy t^ise ~ tdecay — ^d/rg- Those times might be 
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comparable to the light crossing time of the blob if dissipation commences at a radius ~ 
as in the case shown in [J] and [5l In the case of delayed deceleration the rise and decay times 
are determined by the delay: — tdecay — {^dec — ^d)(l ~ Pv)- The lightcurve tends to be 

more symmetric in this case (see figs [5] for a comparison). The change in the emissivity during 
the coasting (pre-deceleration) phase, as seen in the lower right panel in figure 01 is due to the 
dependence of the intensity of target photons and the density of emitting electrons on r [jsc oc 
for the conical flow considered in the above examples, with intensity profile fs = x~‘^ in Equation 
dH)). This dependence may slightly affect the shape of the lightcurve. 

The shape of the lightcurve displayed in the lower left panel of figure 0] is in qualitative agree¬ 
ment with those computed in Joshi et al. (2014), and is quite typical to observed gamma-ray 
flares. The rough symmetry of the light curve is due to the assumed uniformity of the emissivity. 
In cases where the acceleration of the electrons is conhned to a small region inside the shell, e.g., 
acceleration in shock fronts, we anticipate a faster rise and a slower decay, particularly at very high 
energies, at which the cooling time is much shorter than the light crossing time of the shell. 


3. Discussion 

We considered the dynamics of a dissipative shell in the presence of a strong radiative drag. 
Our analysis indicates that for rapid dissipation substantial deceleration is anticipated in blazars 
on sub-parsec scales (and even parsec scales for sufficiently luminous sources), that should give rise 
to rapid, large amplitude variability owing to changes in the beaming factor. It is worth noting that 
even modest changes in the Lorentz factor can lead to large amplitude variations of the emitted flux, 
owing to its sensitive dependence on the Doppler factor. In principle, this mechanism can produce 
flares with durations as short as the dynamical time of the central engine, at a very high efficiency. 
The radiative drag exerted on the thermal electrons alone should lead to rapid deceleration at 
radii r < 4 x IO^^roT 544 cm, where Tq is the bulk Lorentz factor at the onset of dissipation and 
Lg = I 0 ^^Ls 44 erg s“^ is the luminosity of the ambient radiation field intercepted by the flow, 
provided the dissipation is rapid enough to keep the specific enthalpy large (that is, h » I). This 
is, for instance, the situation in internal shocks, as shown in section 12. 1. 21 but may also occur in 
other cases, e.g., effective magnetic field dissipation in Poynting flux dominated jets (Sikora et al. 
1996). Scattering of ambient photons by non-thermal electrons accelerated in situ will dominate 
the radiative force if the proper scale over which electrons are accelerated exceeds ct'^/^e: where 
is the cooling time of the non-thermal electrons, as measured in the rest frame of the flow, and 
is the fraction of total energy carried by the non-thermal population. 

In general, the effect of radiative drag is expected to be more prominent in FSRQs than in 
BL Lacs. Nonetheless, Equation (1191) indicates that deceleration of dissipative shells may also be 
relevant to low luminosity sources, provided that electrons can be effectively accelerated to the 
cooling cutoff. This seems to be case in TeV blazars. It has been shown elsewhere (Levinson 
2007) that if the TeV spectrum emitted during strong flares extends to energies at which the pair 
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production opacity exceeds unity deceleration should be effective. 

It is naively anticipated that the gamma-ray emission produced through IC scattering of am¬ 
bient photons by the nonthermal electrons accelerated in the blob will be correlated with lower 
energy emission generated by synchrotron cooling of the same electrons. However, synchrotron 
self-absorption may give rise to a strong suppression of the radio emission in cases where the de¬ 
celeration occurs well below the radio core. In that case, the gamma ray flare will either precede 
the ejection of a superluminal component, or not be accompanied by one at all, depending on the 
asymptotic bulk Lorentz factor and the synchrotron cooling time. On the other hand, the onset of 
dissipation depends on the duty cycle (the time interval between ejections of consecutive shells in 
the case of internal shocks), and is expected to occur over a range of scales, even in an individual 
object. Blobs that dissipate their energy at large enough radii will not experience strong decelera¬ 
tion. As a consequence, a variety in the behavior of flares is expected, as indeed revealed by recent 
multi-waveband studies (e.g., Marscher et al. 2011). 

We thank Noemie Globus for help. This research was supported by a grant from the Israel 
Science Foundation no. 1277/13. 


A. Derivation of the source terms 


The general expression for the source term associated with Compton drag is (Phinney 1982; 
Sikora et al. 1996, Van Putten & Levinson 2012): 


= -aTj^J ^frfePuk’' 


k^ + 


{p„k’')p^^ 


mi 


(Al) 


Here fr and fe denote the distribution functions of the photons and electrons, respectively. In the 
rest frame of the blob the electron distribution /g is isotropic, and it is convenient to compute the 
source terms there, and then to transform. We have 



—p^k'^ +p' ■ k' 


-p'^k'^ + p'k'p = -p'°A:'°(I - 


(A2) 


where p is the cosine of the angle between the photon and electron directions, and p'^ jme = 7, 
p''^/me = 7/3*. The zeroth component reads: 


f d^p' , 

"V P-" j 

f d^k' ,, ,0 {p'uk'np^ 

JrJePv''' 1 ^2 

(A3) 

OT J 1 

^ dyr.k!^ [(1-/3m)-7'(I-/3^)2] 

(A4) 

CO 

o 

/ j dp'p'^f^ j dp [(I - f3p) - 7^(1 - (3pf] . 

(AS) 



(A6) 
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Now, n'g 


/ d^p'fe 

5 ;° 


2 f dp'p'‘^fe, and using f dp(l — = 2(1 + /3^/3) one has 

-^ar I d^k'k'^f, I = -^arn'g < 7^/32 > 


The spatial component reads; 

ji I {p'uk''^)p'^ 

' ' 9 

= <^tJ d^k'f' j dh'f'ek'^ [(1 - Up) - 7 '/ 3(1 - 

= or y d^k'k'^fl. J dp'p"^f^(2 + ^7^/5^) = cttK < 1 + 27^/3^/3 > 


5'* = -UT 


f d?k' 


P 


,/o 


fco 


frfePlk‘ 


(A7) 


(AS) 

(A9) 

(AlO) 


where 

= J ^k'>^k'’'f' (All) 

Now, in the black hole frame the radiation field is isotropic and we have = Ug, T®* = 0, 
= {us/3)Sij. Then 

r;°* = A°Atr^'^ = -V/3fu,. (A12) 

o 

Thus, we finally have for the source terms in the comoving frame: 

S'c = -^o-TUg < 7^/3^ > Usr2(l + /3r/3), (A13) 

S'J = < 1 + 272/32/3 > Ugr2/3f. (A14) 

o 

To obtain the source term in the Lab frame we perform a Lorentz transformation, 

5° = + ^/3^5^ = -laTn'gr\,[< 72/32 > (1 + /32) + /32]. (A15) 

SI = r/3f 5g° + = -^aTn'gr3/3fu,[l + (5 + /3^) < 72/32/3 >] (A16) 

o 


For cold electrons < 72/32 >= 0, and S'g® = 0, as required, since in the Thomson limit invoked here 
the scattering is fully elastic and there should be no energy loss in the rest frame of the blob if the 
electrons are cold. In the Lab frame we then recover the result S'g = — |cJTUg T^/3^Us. 


B. Flow equations 

In spherical coordinates, the radial expansion of a hydrodynamic shell is governed by the 
equations 


dt{wT^-p) + ^dr{r^wT^(3r) = 5g°, 

(Bl) 

dti'wT'^Pr) + ^drir'^wT'^jS^) + drP = 5/, 

(B2) 
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where 5° and are source terms that account, respectively, for energy and momentum losses by 
the radiative drag. The above equations can be combined to yield 

wr ^^ + Mp + drP = s: - /3r5°, (B3) 

in terms of the Lagrangian derivative d/dt = dt + Prdr- 

We suppose that at time t the dissipative shell is contained between the radii ri (t) and r 2 (t) = 
ri(t) + A(t). For simplicity we suppose that the length of the shell A is constant in the black 
hole frame. The velocity of the shell is then uniform with /3r(i) = dr 2 /dt. Transforming to the 
coordinates T{r,t) = t, C(b^) = ^2(0 — i", we have 

dr = -9c > 

dt = dr + I3rd(;, 

^ = dt + ^rdr = dr- 
dt 

Using the relation /3rd/3r = T^^dT, Equation ()B3I) gives 

to order 0(r“^). Equation ()B1D yields 

— p) — wV'^/dr — dcp — wdr InT = S^. (B8) 

dt r 

To order 0(A/r2) we have 2/3r/r = {2/r2)dr2/dt = r^'^dr'^/dt, and the latter equation reduces to 
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Fig. 1. — Lorentz factor profiles computed numerically using Equations (jl7p - (jl8p with cooling 
limited injection (Eq. [9|), and different values of a, as indicated. 



Eig. 2.— Pair-production optical depth as a function of photon energy . 
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Fig. 3.— Spectral energy distribution of ERC radiation emitted from a decelerating blob following 
onset of dissipation, for a non-delayed deceleration (left panel) and delayed deceleration with r^ec — 
Td = A/(l — /3r) (right panel). 
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Fig. 4.— Sample light curves for non-delayed deceleration (upper panels) and delayed deceleration 
with Vdec ~ = A/(l — /3r) (lower panels). The right panel in each case exhibits the (normalized) 

flux at a photon energy of 1 GeV measured by a distant observer as a function of observer time, 
and the left panel the corresponding emissivity as a function of radius r = ct. The two curves in 
the upper left panel correspond to different jet opening angles, as indicated. 
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Fig. 5.— The effect of delayed deceleration on the shape of the light curve, 
in all cases shown is 6j = 7r/2. 
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